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We have investigated the dependence of galaxy clustering on their stellar mass at z~l, using 
the data from the VIMOS-VLT Deep Survey (VVDS). We have measured the projected two- 
point correlation function of galaxies, Wp{rp) for a set of stellar mass selected samples at an 
effective redshift < z >= 0.85. We have control and quantify all effects on galaxy clustering 
due to the incompleteness of our low mass samples. We find that more massive galaxies are 
more clustered. When compared to similar results at 2; ~ 0.1 in the SDSS, we observed no 
evolution of the projected correlation function for massive galaxies. These objects present a 
stronger linear bias at z ~ 1 with respect to low mass galaxies. As expected, massive objects 
at high redshift are found in the highest pics of the dark matter density field. 



1 Introduction 

In the currently accepted scenario, galaxies are thought to form within extended dark-matter 
haloffi, which grow through subsequent mergers in a hierarchical fashion. A major challenge in 
testing this general picture is to connect the observable properties of galaxies to those of the 
dark-matter halos in which they are embedded, as predicted, e.g., by large n-body simulationiP. 
Recent theoretical works seem to indicate that a fairly direct relationship indeed exists between 
global galaxy properties (e.g. their stellar or total baryoni c mass, o r their luminosity) and the 
halo mass, before it is accreted by a larger dark-matter 

At the current epoch and at high redshift, the observed clustering of galaxies is found 
to depe nd significa ntlv o n their specific properties, such as luminosit>*^ ^ ^ , color or spectral 
typ^^, morphologj^^^ and stellar masP^. The dependence of clustering on stellar mass in 
the local Universe has been shown in the Local Universe, from a large sample of galaxies of the 
Sloan Digital Sky Surve;^^. Galaxies of larger mass are seen to be more clustered than low-mass 
ones, with the effect increasing above the characteristic knee value M* of the Schechter mass 
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function. Using the VIMOS-VLT Deep Survey (VVDS), we have investigated for the first time 
the evolution with redshift of the galaxy clustering as a function of stellar mass at 2 ~ i^. 

The VVDS is performed with the VIMOS multi-object spectrograph at the ESO Very Large 
Telescop( ^^ and corn plemented with multi-color imaging data obtained at the CFHT and NTT 
telescopePZIISIinEQ]. We have used the so-called "first epoch" data, collected in the VVDS-02h 
"Deep" field. This is a purely magnitude limited sample to Iab = 24, covering an area of 0.49 
square degrees with a mean sampling of ~ 23%. The full sample includes ^ 11000 redshifts. 
Details about observations, data reduction, redshift measurement and quality assessment can 



be found elsewher(^lj. 

We use a Concordance Cosmology with Q^n = 0.3 and = 0.7. The Hubble constant is 
normally parametrised via h = //q/IOO to ease comparison with previous works. Stellar masses 
are quoted in unit of /i = 1. All correlation length values are quoted in comoving coordinates. 



2 Stellar mass and incompleteness 



Stellar masses for all galaxies in the VVDS catalogue were estimated by fitting their Spectral 
Energy Distribution, as sampled by the VVDS multi-band photometry, with libraries of stellar 
population modelP^l. The robustness of the method and intrinsic errors have been fully testeePS^. 
The typical error on stellar mass is ~0.1 dex. 

We concentrated our study of the dependence of galaxy clustering on stellar mass in the 
redshift range [0.5-1.2]. Such a sample includes 3218 galaxies more massive than 10^ Mq. 
However, the VVDS is a flux limited survey, and even if the luminosity is correlated to the 
stellar mass, this relation has some scatter and then, some faint galaxies might be missed in 
a given stellar mass range. We investigated the completeness limit using the data themselves 
and 40 mock catalogues complete in stellar mass built from the Millennium run dark-matter 
output^'^ coupled to semi-analytical modelP^. 
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Figure 1: Stellar mass completeness at Iab = 24 as a function of redshift estimated from 40 mock catalogues. The 
green area indicates regions where all galaxies are fainter than the apparent luminosity thresold. The blue area 
shows region where no galaxy was found in the mock at any Iab apparent magnitude. The red points correspond 

to the VVDS data. 



As an illustration, Fig. [T] shows the stellar mass vs. redshift relation of the mock catalogues, 
with, coded in grey scales, the fraction of galaxies brighter than I^b = 24 in narrow redshift 
(Az = 0.01) and stellar mass (Alog(M) = 0.01) bins. According to these numerical simulations, 
we miss ~ 30% of galaxies in the stellar mass range log(M/MQ) = [9.5 — 10] between redshift 
z=0.5 and z=1.2. Samples of galaxies more massive than lO^'^ Mq are almost complete. It is 
interesting to note that, rather realistically, the typical galaxies producing the incompleteness 
at faint fluxes in the mock samples are red objects with high mass-to-light ratio. 



3 Galaxy clustering as a function of stellar mass 



We measure the two-point correlation function ^(rp, vr) and its projection along the line of sight 
Wp{rp^^ using standard estimatoP^. We take into account the complex observing strategy of 
the survey. We checked using the mock catalogues the effect of stellar mass incompleteness 
on the clustering measurement, comparing the projected correlation function for mock samples 
complete in stellar mass and selected at Iab < 24. We concluded that, when incompleteness 
is the strongest (i.e. for log(M/M0) = [9 — 9.5]), the amplitude of the correlation function is 
underestimated on small scales by a factor up to ff^. 

Taking into account these effects in our conclusions, we measured the projected correlation 
function for several stellar mass samples from the real data at z ~ 1. We observed for the first 
time at this redshift that more massive galaxies are more clustered than low mass galaxies, the 
correlation length tq increasing from 2.76~!:o\l h"^Mpc for log(M/M0) > 9 to 4.28^0 45 h'^Mpc 
for log(M/M0) > 10.5. This dependence on stellar mass was predicted b y h ydrodynamical 
simulationPS] and is also observed in the mock catalogues we have been usin^^. 

We then compared our measurement to similar analysis made at 2; ~ 0.1 with SDSS dati^. 
Figure [2] shows a faster evolution of the correlation function Wp{rp) for galaxy less massive than 
"I^qIO.s ]^^^ while the clustering properties of massive galaxies do not appear to evolve with time. 
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Figure 2: Comparison of the projected correlation function Wpifp) from SDSS data at z~0.15 {black cross) to 
VVDS measurement at Z'^O.SS {magenta filled dots and solid lines) in 3 ranges of stellar masses. The green dashed 
line is a power-law reference line drawn with ro — 5 h~'^.Mpc and 7 — 1.8. Error bars on VVDS measurement have 
been estimated from the variance among 40 mock catalogues. Blue open circles indicate VVDS measurements 

corrected for stellar mass incompleteness. 



Thereby, the most massive galaxies display an evolution of their linear bias factor (from 
bi = 1-62 it 0.18 at z~l to bi = 1-40 it 0.03 at z~0.1), while it remains roughly constant for 
lower mass objects (6l ~ 1.3). Massives galaxies are found at z ~ 1 in the highest peaks 



of the density This finding is expected in a hierarchical scenario in which the most 

massive peaks of the mass density field collapse earlier and evolve fasteP^l. This interpretation 
qualitatively supports a scenario in which the stellar mass of a galaxy is essentially proportional 
to the mas s of th e dark- matter halo in which it was last the central object, consistent with recent 
simulationP2II3]. 
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